General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 





i>nr<A 


51 


MSC - (M217 


NATIONAL AEIONAUTICS AND SPACE ADMINISTRATION- 


PROJECT/SPACE SHUTTLE 

SPACE SHUTTLE GUIDANCE, 
NAVIGATION AND CONTROL 
DESIGN EQUATIONS 




VOLUME IV 

DEORBIT AND ATMOSPHERIC 
OPERATIONS 












APRIL 15, 1971 

5 - C^CODE) 

2 O X' 3 ^- 

4 ^ ^CAIEGORV! 

I ,1^',AC.oVwXO»AONUMI>t'<l 


SYSTEMS ANALYSIS BRANCH 
GUIDANCE AND CONTROL DIVISION 

MANNED SPACECRAFT CENTER 

HOUSTON.TEXAS 


^ * -I - S- 1 


MSC - 04217 


NASA SPACE SHUTTLE PROGRAM WORKING PAPER 


SPACE SHUTTLE GUIDANTc, NAVIGATION 

AND CON’ ; OL 
DESIGN EQL.'-TIONS 


VOLUME IV 

DEORBIT AND ATMOSPHERIC OPERATIONS 


April 15, 1971 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
MANNED SPACECRAFT CENTER 
HOUSTON. TEXAS 


Prepared by 

Systems Analysis Branch 
Guidance and Control Division 


k . r, V 

K. J. Cox, Chief 
Systems Analysis Branch 


Authorized . / 

Distribution 

^^'LMaxime A. Foget / 

' Director of Engineering and Development 


% 


TABLE OF CONTENTS 


VOLUME I 

1. PJRFOSE 1-1 

2. SCOPE 2-1 

3. APPLICABILITY 3-1 

4. DEFINITION OF TERMS 4-1 

5. REFERENCE DOCUMENTS 5-1 

6. GN&C SYSTEM DESCRIPTION 6-1 

7. GN&C SOFTWARE FUNCTIONAL REQUIREMENTS 7-1 

8. COORDINATE SYSTEMS AND TRANSFORMATIONS 8-1 

VOLUME II - PREFLIGHT THROUGH ORBIT INSERTION 

9. DESCRIPTION OF EQUATIONS 9-1 

9.1 Preflight 9.1-1 

9.1.1 GN&C System Initialization, Monitor, Test 

and Checkout 9.1-1 

9.1.2 Prelaunch Targeting 9.1-1 

9.1.3 Sensor Calibration and Alignment 9.1-1 

9.2 Boost Phase Program 9.2-1 

9.2.1 Rapid Real-Time State Advancement During 

Specific Force Sensing 9.2-4 

9.2.2 Boost Monitor 9.2-15 

9.3 Separation 9.3-1 

9.3.1 Attitude Control 9.3-1 

9.4 Orbit Insertion Program 9.4-1 

9.4.1 Navigation 9.4-1 

9.4.2 Powered Ascent Guidance 9.4-3 

9.4.3 Attitude Control 9.4-25 

iii 


TABLE OF CONTENTS (CONTINUED) 




PAGE 

9.5 Ascent Aborts 9.5-1 

9.5.1 Abort Decision and Mode Determination .... 9.5-1 

9.5.2 Abort Targeting 9.5-1 

9.5.3 Abort Guidance 9.5-1 

9.5.4 Navigation 9.5-1 

9.5.5 Attitude Control 9.5-1 

VOLUME III - ORBITAL OPERATIONS 

9.6 Orbital Coast 9.6-1 

9.6.1 Navigation 9.6-1 

9. 6. 1.1 Conic State Extrapolation 9.6-3 

9. 6. 1.2 Precision State and Filter Weighting 

Matrix Extrapolation 9.6-33 

9.6.2 Sensor Calibration and Alignment 9.6-62 

9.6.3 Orbital Coast Attitude Control 9.6.62 

9.7 Orbital Powered Flight 9.7-1 

9.7.1 Targeting 9.7-2 

9.7.2 Navigation 9.7-2 

9.7.3 Guidance 9.7-2 

9.7.4 Attitude Control 9.7-2 

9.8 Rendezvous Mission Phase 9.8-1 

9.8.1 Rendezvous Targeting 9.8-3 

9.8.2 Relative State Updating 9.8-90 

9.8.3 Rendezvous Guidance 9.8-132 

9.8.4 Rendezvous Attitude Control 9.8-132 

9.9 Station Keeping Mission Phase . . 9.9-1 

9.9.1 Relative State Estimation 9.9-2 

9.9.2 Station Keeping Guidance 9.9-2 

9.9.3 Station Keeping Attitude Control 9.9-2 

iv 


« 

/ 


TABLE OF CONTENTS (CONTINUED) 

PAGE 

9.10 Ducking and Undocking 9.10-1 

9.10.1 Docking and Undocking Navigation 9.10-5 

9.10.2 Automatic Docking Control Law 9.10-32 

9.11 Docked Operations 9.11-1 

APPENDIX A - APPLICABLE APOLLO SOFTWARE A-1 

VOLUMF IV - DEORBIT AND ATMOSPHERIC OPERATIONS 

9.12 Deorbit and Entry 9.12-1 

9.12.1 Deorbit Targeting 9.12-1 

9.12.2 Entry Navigation 9.12-1 

9.12.3 Entry Guidance 9.12-4 

9.12.4 Reentry Autopilot 9.12-16 

9.13 Transition 9.13-1 

9.13.1 Transition Attitude Control 9.13-2 

9.14 Cruise and Ferry Cruise 9.14-1 

9.14.1 Navigation 9.14-1 

9.14.2 Guidance 9.14-1 

9.14.3 Attitude Control 9.14-1 

9.15 Approach and Landing 9.] 5-1 

9.15.1 Navigation 9.15-1 

9.15.2 Guidance 9.15-1 

9.15.3 Control 9.15-1 

9.16 Horizontal Takeoff 9.16-1 

9. L7 Communications and Pointing 9.16-1 

9.18 Failure Detection 9.16-1 

VOLUME V - FLOW DIAGRAMS 

10. FLOW DIAGRAMS 


v 




TABLE OF CONTENTS (CONTINUED) 


PAGE 

VOLUME VI - CONSTANTS AND KEYBOARD 
ACCESSIBLE PARAMETERS 

11. SUMMARY OF CONSTANTS 

12. REQUIRED KEYBOARD ACCESSIBLE PARAMETERS 


Vi 


% 


PURPOSE 

The purpose of this document is to specify the equations 
necessary to perform the guidance, navigation and control onboard computa- 
tion functions for the space shuttle orblter vehicle. This equations 
document will: (1) establish more specifically, than on a functional 

level, the GN&C computational requirements for computer sizing, (2) provide 
GN&C design equations specification to develop demonstration software for 
hardware feasibility testing, and (3) define the hardware Interface 
requirements with the GN&C subsystem software. The document will provide 
a standard of coiranunlcatlon of Information concerning the GN&C equations, 
and will provide a means of coordination of GN&C equation development. 
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SCOPE 


« 


This is Volume IV of the document which defines the Guidance, 
Navigation and Control (GN&C) design equations sequencing and Interfaces 
for the computations required In the GN&C Subsystem for all mission 
phases of the Shuttle Orblter flight. The equations are Intended to 
satisfy the functional requirements specified In Reference a. This 
document will describe In mathematical, logical, and opiratlonal language 
all the details necessary to Initiate and carry out the design of the 
required computer modules (subprograms) for the GN&C functions. 

The document will be organized Into six volumes. Volume I 
contains Sections 1 through 8, which provide Introductory Information for 
the document. Volume II contains the current detailed equations for the 
preflight, boost, separation, orbit Insertion and ascent abort phases of 
the Orblter operation. Detailed equations for orbital operations of the 
Orblter, which Include the orbital coast, orbital powered flight, rendez- 
vous, station keeping, docking and undocking, and docked operations phases, 
are presented In Volume III. Volume IV contains the current detailed 
equations for the deorblt and entry, transition, cruise and ferry, approach 
and landing, and horizontal takeoff phases of the Orblter. Also, this 
volume will contain the equations for communications and pointing functions 
and the failure detertion function. For this Issue of the document, only 
Volumes 1 through IV are being published. In future Issues, Volume V will 
contain the detailed flow diagrams for the equations. For the Initial 
issue, the flow diagrams for the approved equations are Included with the 
equations In Volumes II, III or IV. For future Issues, the constants used 
in the equations will be summarized In Section 11 and the GN&C parameters 
and variables which can be entered or called via the keyboard will be 
enumerated In Section 12. These two sections will be contained In 
Volume VI of the document. 
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APPLICABILITY 
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This document Is applicable to the Guidance, Navigation, and 
Control (GN&C) Subsystem of the Electronics System of the Space Shuttle 
Orblter Vehicle. It Is applicable to the definition of the shuttle 
computational requirements for the subsystem listed above. It Is applicable 
to the Phase B and Phase C subsystem development. It defines the Manned 
Spacecraft Center Guidance and Control Division Inhouse study baseline 
equations design. 
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9. DESCRIPTIONS OF EQUATIONS 

Thfc detailed equations for the GN&C functions are defined in this 
section. The organization of this section is tentative and will be 
modified so as to present the equations as they are designed in as clear 
a fashion as possible. As an introduction to each major subsection 
(usually a mission phase) , the general CN&C software functions to be 
implemented will be identified and, where appropriate, a conceptual 
discussion and top level flow of the computations, inputs .ind outputs 
will be included in order to understand and summarize what is to be 
covered. This should be an order of magnitude less detailed than the 
flow diagrams of the equations which come later. 

A GN&C Equation Submittal sheet will Introduce each of the GN&C 
equation submittals and summarize the GN&C functions, and identify the 
source and NASA contact for each. 

The detailed data to be presented for each GN&C function within each 
of the major subsections (usually a mission phase) is summarized below. 
Although items 6 through 10 are to be referenced only in the equations 
document, they are required submittals before the equations can be 
approved and finalized for flight software development. 

1. Functional Requirements 

The specific functional requirements (from the GN&C 
Software Functional Requirements Document) which are 
satisfied by the equations should be identified. 

2 . Functional Diagram 

A brief functional explanation and description of the 
overall concept and approach. A functional block 
diagram should be used where clarity is enhanced. 

Inputs, outputs, and interfaces will be provided. 

3 . Equations and Flows 

Detailed equations and a descriptive text which guides 
the reader through the flows of Section 10 should be 
provided. The minimum frequency of the computations 
shall be specified and rationale given or referenced. 
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U . Coord Irinle System 

Thf coordinate syatcms used shall be defined. 

5 . Co nsta nt s/Voriablcs Su mmary 

Constants and variables shall be summarized In tabular 
form with the following Information: 

a. Variables/constants symbols and definitions 

b. Units 

c. Allowable quantization 

d. Range of values 

6. FO RTRAN Coding 

The FORTRAN coding of the function for verification using 
the Space Shuttle Flight Simulation (SSFS) will be 
referenced. 

7. Simulation 


The SSFS r.pecif Icatlons , description and user's guide 
used to verify each GN&C function will be referenced. 

8 . Testing 

Test plans and test results will be referenced. 

9 . Derivatio n 

The mathematical derivation of the equations including 
all mathematical assumptions shall be referenced. 

10. Assumptions 

The following will be referenced: 

a. Avionics baseline system assumed 

b. Reference missions assumed 

c. Vehicle mass properties assumed 

d. Propulsion models assumed 

e. Environment models assumed 

f. Error models assumed 

The major subsections of this section are identified and partially 
expanded in the following. 
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9.12 OrorBIT AND ENTRY 

The Deorblt and Entry phase of the Shuttle mission begins with 
the Shuttle pertorr.ing the Deorblt Targeting for the Deorblt (or Deorblt 
phasing, if required) burn and terminates at the beginning of transition 
to low angle of attack. The primary landing point and the required orbit 
maneuvers have been determined from pre-phase planning computations. The 
GN&C Software functions to be performed during this phase are; 

1. Deorblt Targeting will be performed in the coast 
period prior to the deorbit burn(s). Targeting 
outputs (for a specified landing site input) will 
Include parameters as tlme(s) of ignition, Initial 
attitude for the burn(s), burn time, total AV 

for each burn, IMU realignment schedule, time of 
reentry, transition time, cross and downrange re- 
quired, Jet engines ignition time, and time of 
touchdown. 

2. Orbital navigation using ground beacon or star/ 
horizon measurements as required to update Shuttle 
State prior to thrusting maneuvers and aerodynamic 
entry augmented by available external measurements. 

3. Powered Flight guidance resulting in commands to 
autopilot (attitude), thrust vector controls, or 
thrust throttle commands to achieve desired con- 
ditions (state and attitude) for transition and 
limit entry g-loadlng and heating. 

4. Autopilot computations to compute reaction control 
system commands to achieve the proper attitude. 

Figure 1 displays a functional flow diagram of the Deorbit and Entry 
Software Functions. 


9.12.1 DEORBIT TARGETING (TBD) 

9.12.2 ENTRY NAVIGATION (TBD) 
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GN&C SOFTWARE EQ UATION SUBMITTAL 


Software Equation Section: Reentry Guidance Submittal No. 2 

Function; Provide guidance commands during entry 

Module No. 0G5 Function No. 9, 10 >11 (MSC 03690) 

Submitted by: John W. Sunkel Co. GCD 

Date : Dec 1970 

NASA Contact: W. Peters Organization; GCD 

Approved by Panel III: |^ . 3" . CvfC Date: / — 

Summary Description; A fast-tlme-predlctlon (FTP) guidance scheme Is used 
which Is based on sequential rapid Integrations of the equations of motion . 
Senaitivities of all the future states to a control change are known from 
perturbation solutions. The predi ctor Is a fourth-order Adams -Bashforth — 
formula: the correction Is a fourth-order Ad ams Moulton formula. 

Shuttle Configuration: Not provided; however, the general concept should 

be Independent of shuttle configuration. 

Comments : 

(Design Status) Coded and In test. 

(Verification Status) Preliminary test results provided. 

Panel Comments: 
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1. 


Introduction 


The Fast-Time Prediction (FTP) Guidance method documented In 
refs. (1) and (2) applied only to flight In a vertical plane over a 
spherical nonrotating earth. The method has now been upgraded to Include 
a more accurate dynamic model as well as predictive crossrange logic. The 
significant revisions are described In the next three sections of this memo. 
Reference 4 presents seme typical results obtained during checkout of 
the new routine. Table 1 defines nomenclature used in the memorandum. 


2 . Equations of Motion 


FTP guidance Is based on sequential rapid Integrations of the 
differential equations of motion-i The previous version of the FTP 
guidance routine used Chapman's equation which applies only to flight In 
a vertical plane and in regions where drag is significant and where the 
flightpath angle is small. Furthermore, coriolis terms are not Included. 
These restrictions have been removed in the FTP-4D program without 
significantly Increasing the complexity of the equations that are 
Integrated. Velocity Is used as the Independent variable to eliminate 
one equation. In addition, the set Is nondlmenslonallzed since. In 
general, well-conditioned equations can be Integrated with larger step- 
sizes. Therefore, the equations that are integrated in the FTP-4D 
guidance program take the following forms: 
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Table 1 


DEFINITION OF TERMS 

Dimensions 

Variable (ND*nondimens ional ) 


6 

scale height in exponential 
atmosphere 

ff^ 

BIAS 

The control parameter which 
increases or decreases the 
in- )lanc L/D 

ND 

Cd 

dra.i coefficient 

ND 

f(?) 

a prespecified function deter- 
mining the shape of the in-plane 
component of L/D 

ND 

A 

g 

gra /itat ional acceleration 

ND 

A 

Y 

flight path angle 

ND 

f; 

alt Itude 

ND 

• 

1 

instantaneous trajectory inclina- 
tion angle 

ND 

A 

i 

crossrange distance (♦ to right 
of track) 

ND 

L/D 

maximum lift-to-drag ratio 

ND 

M. 

non di mens i onal i zing -mass 

4.21x10^ si 

y 

angle of latitude 

ND 

A 

angle bank about velocity vector 
(”♦• rolled right) 

ND 

A 

azimuth with respect to downrange 
direction (•*• to right of track) 

ND 

• 

Q 

heating rate 

Btu/f t ^ /sec 

A 

r 

distance from earth center 

ND 
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Variable 

Rq nondimensionalizing length 
earth radius 


assumed sea-level atmospheric 
density in exponential 
atmosphere . 

reference area 

V velocity relative to a rotating 

atmosphere 

Vq nondimensionalizing speed 

earth spin rate 
e * 

X downrange distance 


Dimens i ons 
(ND*nondiinensional) 

a.issxio"^ ft. 

ND 

ND 

ND 

ND 

2.57x10^ ft/se.. 

ND 

ND 
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FACTORIZATION TERMS 


T1 


I 


T2 

T3 

T5 

T6 

T7 

T8 

T9 

TIO 

Til 

T12 

T13 

T14 

T15 

T17 

T19 

T20 

T21 

T22 

T23 

T24 


u) cos i 
e 

-6hR, 

e 

cos Y 
sin 
cos $ 
sin } 
cos J 
sin ^ 

(T24)/(T19) 

^ 2 

- (in) (TS) (T30) 

- (Til) [ (Tl) (T3) (T15) (T12) + (T10)(T21)] 

(B)(R(,) 

g/v 

- (Tl) (T3) (T12) - (TIO) (T30) 
r ♦ h 

e 

(T2)^ 

r/v 

1/CT22) 

1/(T19) 
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T25 

T26 

T27 

T28 

T29 

T30 


f^(T23) 

T17 - (T22)(T21) 
(T1)(T3)(C)(T7) 
(TIO) (T25) 

(T25) (T5) 
g - (T20)(T21) 
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jj sin 4) 
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sin i sin 


(4) 

A 

dl 

— Sr 

dv 

^5 

V ' 

jV sin 
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Factorization of these equations using the terms defined In the appendix 
results In the following equivalent equations: 


“ (T24) Rt27) - (T5)(T26) + 2(T2)J 


db 

dv 


(T24) (V) (TIO) 


dx 

57 


(T24) (T25) (T5) 


dv 


(T24) (1'6) jj + 2 sin i sin J 


dl 

dv 


(T24) fv)(T8) 


(la) 

(2a) 

(3a) 

(4a) 

(5a) 


The bank angle 4),used in the integration of equations (1) - (3) Is obtained 
from .the exoression 
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-1 

COS 


f('?0 BIAS 

l7T) 


( 6 ) 


where f (v) is a predetermined function that specifies the shape of the 

A 

vertical plane component of L/D as a function of v. BIAS is the vertical 

A 

shift of the L/D cos ^ curve obtained from the previous guidance cycle. 
Thus each new guidance command is followed by a prediction of the range 
error using that command. 
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9.12.3 Entry Guidance (cont’d) 

The gain used in the dovnrange guidance Is obtained by Integrating the 
set of perturbation equations associated with equations (1) - (3). These 
equations are: 


d 3 y 

dv 3 (L/D cos ^ 


d 3h 

dv 3(L/D cos (f) 


d^ 3x 

^ 3(L/D cos <p) 


3x 3(L/D cos 4^) 


oil 

b cos Q J 

X 


3(IVn cos J) 


sfi 


(L/I) cos 



3(L/1-' cos ^3 


As 

dh 

3fi 

3(L/I^ cos 4) 


* / • " 


Three partial derivatives Indicated in equations (7) - (9) can be factored 
and simplified considerably before programming to reduce computation time. 
The resulting forms are: 

- (T2l.) [(-TS)(T13) + (T24)rn0)] f (T1 3 H 7) 1./ D 

p Y 

^ 2(T2}] (1 


{(T27)L/I) - (TS)(T2C>) 2(T2)}(T1.;3 

1 

T24 jCT27) (T1S)L/D ^ (T5)(T2])l/v - 


(T2 3) 

U 2 u j ' 


v{(-T24)(T5) (T10)(T13)} 
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3f 

—~K- 

3h 


nTlO) (T14J 


3f, 

3x 
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3f 


(T28) 


TrfuT ■ 


3f 


(T29HT14J - (T24)(T29) 


31' 


3x 
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3f 


3(I./D cos (J)) 


(T24) (T27)/ ,^17) 


(H) 

(1.M 

(H)) 

(17) 

(H)) 


The T-factors in equations (10-(19) are defined in Table 

The gain for the downrange guidance is 1 / ( 

/ \3(L/n cos C 


at 


cutoff altitude. Since equations (7)-(9) are integrated simultaneously 
with equations (l)-(5) on each guidance call, the gain used is always 
based on the latest control history. The result is a very stable 
guidance response exhibiting little or no overshott. 


Integration of equations (l)-(3) and (7)-(9) provides an abun- 
dance of information that is useful. At every guidance call, all future 
state variables are predicted. Thus accurate predictions of the magnitudes 
and locations of the maximum load factor and heating rate are obtained. 
These could be displayed to the crew or used to alter the flightpath if 
their magnitudes are excessive. Note also that the sensitivities of all 
the future states to a control change are known from the perturbation 
solutions. Therefore, if it is necessary to alter the command to avoid 
exceeding a load factor or a heating rate limit, the proper gain is known 
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9.12.3 Kntrv Guidance (cont'd) 

for that control loop. For example, if the maximum heat rate predicted 

. 2 

is Q and the allowable limit if 100 bta«ft jec, the first-order 
max 

estimate of the change in required to avoid the excess could be obtained 
from the expression 


A(L/1> cos 5) « ^’(Qnax ' 

where 


1 

G 


Jx . L9 

Dy 9(L/n cos oj a( 17 I) cci’ <>) 


( 21 ) 


Both 


• • 





cos 

to equations 


■*' — and 

(7)-(9). 


are simple algebraic expressions; 


— are available from the solutions 

8(h/r) cos C) 


3. Integration Method 

To reduce the number of function evaluations, a predictor- 
corrector integration method is used in the FPT-4D guidance routine. 

The predictor is the fourth-order Adams-Bashforth formula and the 
corrector is the fourth-order Adams-Moul ton formula (see ref. 3). Four 
steps of a fourth-order Runge-Kutta method are used to start the method. 
The integration formulae are built into the guidance routine to eliminate 
calls to external modules. The predictor-corrector formulae together 
with nondimensionallzation and proper factorization of the equations 
allows the integration to proceed very rapidly. A preliminary checkout 
run made 311 guidance calls in 111 seconds on the 1108. Simulations to 
date assume that the guidance is called every’ 3 seconds of real time. It 
is probable that lower-order integration formulae and larger integration 
stepsizes could be tolerated if necessary. 

4 . Crossrange Logic 

Currently, integrated crossrange predictions are used as the 
basis for the crossrange guidance logic in the FTP-4D routine. The 
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differential equations used to calculate the crossrange are 


d’^ 

dv 


~ 1 4 P.ve ^ C^S k 


2 ^0 


1) 


+ 2w sin i sin u 
‘ ) 


d 

dv 


1 \ ^ 7 I 

- V SI n i;/ 

• I I 


( 22 ) 
(2 3 ) 


During each guidance pass, equations (22) and (23) are integrated using 
the most recent history multiplied by minus one. Thus the result of 
the integration is a prediction of the crossrange that would result if the 
sign of the roll command were reversed. A five-degree per second reveisfcl 
rate (over the top) is assumed to Improve the accuracy of the predicted 
reverse crossrange capability. The ratio of the vehicle's current 
crossrange to its reverse crossrange capability is used as the basis for 
commanding a roll reversal. The threshold values Increase linearly as 

A A 

V decreases. Typically, the threshold is set at 80% for v = ^8 and 
100% for V =.15. Accurate crossrange capability predictions obtained 
from Integrating equations (22) and (23) allow large percentages of the 
crossrange capability to be 'used up' during each roll segment of the 
trajectory. The result is a minimization of the numbir of roll reversals 
required. Furthermore, the logic does not tend to command reversals in 
rapid succession near the end of the trajectory. 

Several refinements are being examined. Most are directed 
toward Improving the crossrange accuracy. Among these are compensation 
for the lag introduced by finite computation times and improved design 
of the switch curve. As suggested in Ref. (2), a stable downrange channel 
is not seriously affected by reversal commands from the crossrange logic. 
However, the downrange logic strongly affects crossrange accuracies 
because as downrange errors are nulled, crossrange capabilities are 
altered significantly. Several possibilities are being examined to 
minimize the downrange-to-crossrange interference. 
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Software Equation Section D eorblt and Entry Submittal No. 

F unc t Ion Reentry Autopilot 

Module No. QC5 Function No. -1, ~ 3 |[MSC 03690) 

Submitted By; G. L. Zacharlas Co. MSC/GCD 

(Name) 

Date; December 2, 1970 

NASA Contact; W. H. Peters Organizatio n EG2 

(Name) 

Approved by Panel III K. J. Cox D at e Decembei 2, 1970 

(Chairman) 

Summary Description; The objective of the autopilot is to provide 
attitude control of the vehicle from deorbit to transition. Control 
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Reentry Autopilot 


1. INTRODUCTION 

The objective of the Reentry Autopilot is to provide attitude 
control of the vehicle from deorbit to transition. Control in this case 
implies following the roll commands from guidance about the wind velocity 
vector, maintaining the trim angle of attack and maintaining zero side- 
slip angle. No direct control is exercised over the pitch angle, however 
pitch rate damping is used. Roll is controlled by yawing the vehicle. The 
present autopilot has been designed for velocities only as low as Mach 5 
but efforts are being made to extend this to include velocities as low as 
Mach 2, 


Nomenclature 


B 

Ixx, lyy, I-c? 
P. q. r 
(R 
R 
Rc 

Ux , Uy , Uz 
a 


6 

A 

"^ 1 * ^2 
X y 


Scaled value of sideslip angle (B) 

Moments of inertia about principal axes 

Body angular rates about the x, y, z axes 

Roll error 

Roll rate (about 

Roll commands from guidance 

RCS control accelerations (Ux * t /Ixx) 

X 

Relative wind velocity 
Body axis system 
Ajigle of attack 
Trim angle of attack 

Sideslip angle 
Pitch rate deadband 
Computed switch times 

RCS roll torques about body axes 
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Reentry Autopilot (cont*d) 


% 


2. FUNCTIONAL FLOW DIAGRAM 

The chronological sequencing of functions performed by the 
Reentry Autopilot is described in this section and illustrated by the 
functional flow diagram in Figure 1. 

The Reentry Autopilot assumes control after deorbit has been 
ferformed. It is assumed that the trim angle of attack has been achieved. 
The first function performed is to process any roil commands and monitor 
the vehicle roll rate. The necessary roll jet on-times are computed on 
the basis of these inputs. Next, the roll coordinator portion of the 
autopilot is executed. Sideslip control is coordinated with roll control 
and the body roll acceleration U is chosen to aid in a coordinated roll. 

X 

During roll rate hold (i.e., U “0), sideslip and sideslip rate are 

z 

directly controlled in the next portion of the autopilot. Finally, the 
pitch race control commands are computed in the last portion of the 
autopilot . 


3. PROGRAM INPUT-OUTPUT 

The required inputs to the Reentry Autopilot are roll commands 
and actual roll angle, sideslip angle and sideslip angle rate, body 
angular rates, and the RCS control accelerations. It is possible that 
the autopilot may be able to estimate its own accelerations based on 
some parameter such as mass. The outputs of this autopilot are predicted 
jet on-time histories for the various RCS jets. 

Input Parameters 
Body angular pitch rate 
Roll angle about ^ 

Roll angle command 
^dy angular roll rate 
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Reentry Autopilot (cont'd) 
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e 

A 


U^(t) 

Uy(t) 

u,(t) 


RCS control acceleration about Z>body axis 
Sideslip angle 
Sideslip angle rate 
Pitch rate deadband 

Output Parameters 
Predicted X-axls jet firings 
Predicted Y-axls Jet firings 
Predicted Z-axls Jet firings 


4. DESCRIPTION OF EQUATIONS 

The Reentry Autopilot Is divided Into four sections, the roll 
controller, the roll coordinator, sideslip controller and the pitch rate 
damper. The following four sections describe the equations for each of 
the sections of the autopilot. Figure 2 contains a description of the 
coordinate system and variable definition. The system block diagram Is 
Illustrated In Figure 3. 


U P 



f 
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4. 1 Roll Controller 

The following four assumptions are made for the roll dynamics: 

1) The sideslip angle (6) Is zero 

2) The vehicle Is trimmed In angle of attack (a * a^) 

3) There Is negligible Inertia cross-coupling (I « 0) 

X z 

A) There are negligible aerodynamic torques about the 
Z-axls 


Based on the above assumptions and the coordinate system defined 
In Figure 2 , It can be shown that the dynamics are given by: 



Figure 3. System Block Diagram 
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9,12.4 Reentry Autopilot (contM) 

R sin - r 
r - U 

z 

The roll error, Is defined by: 

» ■ (Rc - R) sin 

where Rc is commanded by guidance. 

Substituting the expression for the roll error into the 
dynamic equations yields: 

CR " -r 

r - U 

z 

The above results indicate that an appropriate model for the 
vehicle dynamics is that of a double Integrator. With double integrator 
dynamics, parabolic switch logic may then be used in the (<R, r) phase 
plane. That is, the state trajectory under the Influence of a constant 
acceleration will be a parabola in the phase plane. A typical phase plane 
and trajectory are shown in Figure 4. 

Given an initial error and error rate ((R^, , the phase plane 

switch points may be computed. Then, using the simplified dynamic model- 
ing, the switch times (T^^, T^, T^) may be computed. Thus, the roll 

controller predicts the necessary control history U (t) to drive the 

z 

errors to zero. 
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Figure A. Roll Control Phase Plane 


4 . 2 Roll Coordinator 

A typical phase plane trajectory for a roll maneuver was shown 
In Figure 4. In this example, the Z-axls Jets fired In a positive sense 
from 0 to were off from to T 2 and fired In a negative sense from 
T 2 to T^. This typical thrust-coast-thrust sequence during a roll 
maneuver Is depicted In Figure 5. 



Figure 5. Typical Jet History for Roll Maneuver 
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9.12.4 Reentry Autopilot (cont'd ) 
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Sideslip control Is achieved through the use of a dual-mode 
controller. The two modes of operation are: 

1) When U^^O, the body roll acceleration is chosen in a 
way such as to aid in a coordinated roll. 

2) When U_ ■ 0 (roll rate hold), sideslip and sideslip rate 
are directly controlled. 

The second of these modes is covered in Section 4.3 and the 
first mode is discussed herein. 

A necessary condition for zero sideslip is: 

p(t) - (cot a^) r(t) for all t. 

During the time interval [0, this relation is approximated by choosing 

the pulse width, At, to satisfy: 

p(T^) - p(0) - At 

where p(Tj^) is computed from 

p(T.) - r (T. ) (cot a ) - r (cot a ) 
ill max T 

r is the maximum value of the Z-body angular rate obtained from the 
max 

roll firing. The pulse is then centered about T^/2 as shown in 
Figure o, A similar choice is made for the interval [T 2 , T^]. 



Figure 6. Jet History for Coordinated Roll 
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4 . 3 Sideslip Controller 

During the periods of roll rate hold ('J . 0), sideslip and 

J 

sideslip rate may be directly controlled. A simplified oscillator model 
of the sideslip dynamics may be derived from the linearized lateral 
equations of motion. Restricting attention to the hypersonic regime 
(M > 5), all aerodynamic torques except C ^ are assumed to be zero, 
since those that do not contribute to damping are negligible. Further, 

If the cross product of Inertia neglected, the following approxima- 

tion results: 

2 

6 + u) 6 - U 

where 

5 - (q sin 

U i sin 

By defining the variable ‘B, where 

B = u) B 

the (B, B) phase plane ti ijectorles are circles, and the switch curves 
may be obtained using conventional optimal control techniques. The cost 
function used In determlng the switch curves Is a weighted function of 
the time-optimal and fuel-optimal solutions. The exact shapes and place- 
ment of the curves are functions of the relative weights placed on the 
time and fuel-optimal solutions. 

Figure 7 shows a straight line approximation to a typical optimal 
switch logic for sideslip control; also shown Is a typical phase plane 
trajectory. (Note tht t U = sin Qij) • 
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Given an initial error and error rate 6^), the phase plane 

switch points may be found. Then, using the oscillator modeling, the 
switch times (t^, •••) may be computed. Thus, the sideslip controller 

predicts the necessary control history U (t) to drive (6, 6) to (0, 0). 

A 

4.4 Pitch Rate 

The pitch rate damper control logic is based on the assumption 
that the pitch rate, q, is uncoupled from the lateral dynamics. The 
control law is; 


1) 

If 

|q| i A 

fire Uy to oppose the rotation 

2) 

If 

|q| ‘ A 

do not fire U 

y 




where A is the pitch rate deadband 


Graphical' y, this control law is depicted in Figure 8. Note that the 
pulse time may be predicted from 

At - |q I /Uy 
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9.12.4 Reentry Autopilot (cont’d) 
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Figure 8. Pitch Rate Damper 


5. DETAILED FLOW DIAGRAMS 

This section contains the detailed flow diagrams for the roll 
controller, roll coordinator, sideslip controller and the pitch rate 
damper portions of the Reentry Autopilot. The overall flow diagram 
contained In Figure 1 Indicates how the Individual controllers fit Into 
the overall flow. 


9.12-26 



'A/Ui 


^ Vo/k 




TJ s ^t* 

”^*. " Ti -♦• Atji 
~3 * Ta V Ats 


I Ux(0 = lUx[uJO -Ujl-t)-Ujl*TaVUJt- 


EX IT 


Figure 9. Detailed Flow Diagram, Roll Control 
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9. 12. A Reentry Autopilot (cont*d) 



Figure 10. Detailed Flow Diagram, Roll Coordinator 
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Figure 11a. Detailed Flow Diagram, Sideslip Control 
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9.12.4 Reentry Autopilot (cont*d] 



Figure 11b. Detailed Flow Diagram, Sideslip Control 
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Figure 12. Detailed Flow Diagram, Pitch Rate Damper 
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9.12.4 Reentry Autopilot (cont’d) 

REFERENCES 

1. Zacharlas, G. L., Unpublished presentation "Preliminary Baseline 
Design for Reentry Autopilot." 
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9.13 


TRANSITION 




The transition phase starts at the end of the entry phase Just 
prior to the pltchover from a high angle of attack to tne cruise angle 
of attack. This period of pltchover constitutes the transition phase. 

The transition phase ends with the attainment of the nominal cruise angle 
of attack. Two distinct functions which are performed during this are: 

1. Transition from entry control which is basically 
RCS attitude control to cruising flight control, 
which is basically full aerodynamic control. 

2. Transition from entry angle of attack (high a) 
to cruising angle of attack (low a). 

The GN&C software functions to be performed during this phase 

are : 

1. State advancement during specific force sensing. 

2. Control transition maneuver from high angle of 
attack entry profile to low angle of attack 
cruise profile. 

3. Provide aerodynamic control (RCS backup). 

4. Initiate cruise engine commands. 

3. Initiate cruise navigation and control. 

A flow of software functions during transition appears in 

Figure 1. 


Open Items 

The following items constitute areas where requirements and/or 
systems implementation are unclear and may have significant Impact on 
the GN&C software: 

1. Vehicle stability characteristics during transi- 
tion must be defined. 

2. The control effectiveness of the RCS Jets and the 
aerodynamic surfaces under the conditions en- 
countered from entry to cruise must be defined in 
order to determine their use during transition. 
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9.13 TRANSITION (cont’d) 
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3. The transition sensor configuration must be 
defined. 


Transition Equations 

Baseline equations have not been approved by the GN&C Software 
Equation Formulation and Implementation Panel. The following reference 
is provided to identify the work which has been done. 

1. MIT "Taijk 11 - Digital Autopilot For SSV Orbiter 

Transition," R. F. Stengel, dated December 16, 1970. 

9.13,1 Transition Attitude Control (TBD) 
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9.1A CRUISE AND FERRY CRUISE 

The Cruise and Ferry Cruise phases of the Shuttle apply to the 
atmospheric powered, constant altitude flight conditions prior to final 
descent to terminal landing. This phase includes low altitude earth- 
relative navigation and guidance and aerosurface control. The SW functions 
required in this mission phase are the following: 

1. Powered flight navigation augmented by external 
data to maintain earth-relative position. 

2. Powered flight guidance resulting in thrust level 
control and autopilot (attitude) commands required 
to achieve desired Intersection with tennlnal 
approach navigation for landing. 

3. Autopilot computations to result in proper aero- 
dynamic surface control and trim settings, to 
maintain desired heading and altitude. 

Figure 1 displays a functional flow diagram of the Cruise and Cruise 
Ferry Software Fn ctlons. 


9.14.1 Navigation (TBD) 

9.14.2 Guidance (TBD) 

9.14.3 Attitude Control (TBD) 


9.14-1 


thrust level 
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9.15 


APPROACH AND LANDING 


"Approach" bej^dns as the shuttle advances to the vicinity of 
the landing site with the appropriate position, velocity and attitude for 
landing and has gained acquisition of the ground transponders located at 
the landi.ig site. Landing begins at touchdown and continues until the 
craft reaches ground taxi speed. The GN&C software functions to be per- 
for.:'ed during this phase are; 

1. Approach navigation using measurements from the 

IMU and the ground transponders. Both the shuttle 
state and the landing site state may be estimated. 

?. Approach and landing guidance resulting in commands 

to the autopilot (attitude), thrust throttle commands, 
braking and/or lift flap commands. 

3. Autopilot computations to compute control surface 
commands in order to achieve the command attitude. 

Figure 1 displays a functional flow diagram of the Approach and Landing 
GN&C Software functions. 


Open Items 

The following items constitute areas where requirements and/or 
system implementation are unclear and may have significant Impact on the 
GN&C software: 


1. Approach navigation accuracy requirements have not 
been established. 

2. A ca,' ability to automatically guide during taxi under 
"blind" conditions must be established as a require- 
ment or not a requirement. 


9.15.1 Navigation (TBD) 

9.15.2 Guidanc e (TBD) 

i 

9.15.3 Attitude Control (TBD) 
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HORIZONTAL TAKEOFF (TBD) 


9.17 

9.18 
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COMMUNICATIONS AND POINTING (TBD) 
FAILURE DETECTION (TBD) 
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